Statistics of Fourier Modes of Velocity and Vorticity in Turbulent Flows 
Intermittency and Long-Range Correlations 
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We perform a statistical analysis of experimental fully developed turbulence longitudinal velocity 
data in the Fourier space. We address the controversial issue of statistical intermittency of spatial 
Fourier modes by acting on the finite spectral resolution. We derive a link between velocity structure 
functions and the flatness of Fourier modes thanks to cascade models. Similar statistical behaviors 
are recovered in the analysis of spatial Fourier modes of vorticity obtained in an acoustic scattering 
experiment. We conclude that vorticity is long-range correlated and found more intermittent than 
longitudinal velocity. 
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In fully developed turbulence, most of the experimen- 
tal, numerical and theoretical works Q focus on the 
statistics of the longitudinal velocity increments S r u(x) = 
u(x + r) — u(x). It is now well established that struc- 
ture functions M q (r) = {(5 r u) q ) x , behave as power laws, 



i.e. 



M q (r) 



The universal non-linear evolution 



of C 9 with respect to q is refered to the so-called inter- 
mittency phenomenon: the probability density function 
(PDF) of velocity u is close to Gaussian, while the PDF 
of longitudinal velocity gradients d x u exhibits extremely 
large tails. Another striking property of turbulence is the 
long-range correlation of dissipation events, i.e. (d x ii) , 
up to the velocity correlation length scale L. Many sys- 
tems share same types of behaviors, as financial volatility 
0(a) and electrical transport in granular media 0(b). 
A major issue in turbulence is to derive a possible link 
between long range correlations and vorticity filaments 
0. In the Fourier space, which is an alternative way to 
study turbulence statistics Q , one could expect that the 
Fourier modes of velocity u(k,t) should analogously fol- 
low the same types of behaviors, i.e. (\u(k,t)\ q ) t ~ A: - ^ 9 , 
since k can be interpreted as the inverse of a scale, as 
prescribed in shell models 0. Furthermore, a statistical 
model based on the rapid distortions of the small scales 
predicts strong intermittency in the fc-space 0. Surpris- 
ingly, it is not the case on experimental and numerical 
velocity profiles^as pointed out by the seminal paper of 
Brun & Pumir since the PDF of Fourier modes are 
found to be undistinguishable from Gaussians, whatever 
the value of k. 

However, one of the crucial parameters of the Fourier 
Transform is the finite spectral resolution associated with 
the finite length of velocity profiles hereafter noted £. The 
goal of this letter is to show that, whereas statistics of 
Fourier modes do not depend on k, they depends sig- 
nificantly on the ratio t/L. Firstly, we present such a 
"short time" Fourier analysis of experimental longitudi- 
nal velocity data. Secondly, the statistical dependence 
on t/L is clarified in the context of various turbulent 
cascade models. Finally, we perform a similar analysis 



on experimental data, obtained in an acoustic scattering 
experiment, allowing the direct probing in time of spa- 
tial Fourier modes of vorticity in a turbulent air jet. By 
comparison with the longitudinal velocity data analysis, 
we conclude that vorticity is also long-range correlated 
but more intermittent. 

Let us introduce the Short-Time Fourier analysis 
ue(k,x), which depends on the space variable x and 
the wave vector k and corresponds to the computa- 
tion of the Fourier transform of the longitudinal ve- 
locity u(x) in a window ht of size t, i.e. U((k, x) = 
u(x) ® (e~ lkx ht (x)) , where ® stands for the convolu- 
tion product. We represent in Fig. IHa) the flatness 
of \ui(k,x)\, i.e. Ft(k) = (\u e (k, x)\ 4 )/ (\u e (k, x)| 2 ) 2 , as a 
function of ln(k/k v ), where k v = (tjk) -1 (j\k is the dis- 
sipative Kolmogorov length scale), for several windowing 
length t (see the caption), and for k > 2ir/£. The ex- 
perimental longitudinal velocity signal used in this study 
has been recorded in the Modane's wind tunnel facility 
at a Taylor microscale Reynolds number 1Z\ « 2500 
and thus exhibits a large inertial range. Here, the spec- 
tral window is the Hanning function: hg(x) = cos 2 (7nc/£) 
for x G [—£/2,£/2], he(x) — instead. For inertial wave 
vectors (k < k v ), Ti(k) slightly depends on k but dras- 
tically on £: when £ is of order of several correlation 
lengths {£ 3> L), Ft{k) is close to the Gaussian value 2 
(as found in and when t/L — ► 0, we note a rapid 
increase of the value of this inertial "plateau" . The evo- 
lution of !Fi(k) with respect to t/L is displayed on Fig. 
[Qc) . The error bars have been obtained by a least-square 
fit of the plateaus of the flatness over wave vectors in the 
inertial range. We thus observe that the flatness of the 
Fourier modes behaves as a power law of the scale £, i.e. 
T t {k) - £ a with a = -0.1±0.02. Here, the error bar 0.02 
is large because of the lack of statistics and the fact that 
longitudinal velocity profiles, obtained under the Taylor 
Hypothesis 0, are sensitive to the temporal decorrela- 
tion 0. We have also displayed on Fig. Qfc) the evolu- 
tion of the flatness for an another longitudinal velocity 
profile (i?A ~ 300) obtained in the air jet we will present 
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FIG. 1: (a) Flatness of \ut(k,x)\ as a function of \n(k/k v ) 
(kr, — for different window lengths £ (from top to 

bottom, i/L = 0.057; 0.067; 0.078; 0.089; 0.114; 0.156; 0.200; 
0.228; 0.334; 0.456; 0.91; 1.82; 3.64; 7.28) for the Modane 
velocity signal (1Z\ ~ 2500). (b) Flatness of synthetic velocity 
profiles \u T e ws (k,x)\ (dotted line) and |5? ws (fc, x)\ (solid line) 
vs. ln(fc/fcmax), where fcmax is the Nyquist wave vector and 
for several scales (from top to bottom \og 2 (£ p / L) — -8, -7, -6, 
-5, -4, -3, -2, -1, 0, 1 and 2). (c) Values of the plateaus of 
Figs, (a) and (b) vs. \n(£/L) : (o) Modane velocity data, (□) 
air-jet velocity data (£/L = 0.32; 0.64; 1.28; 2.56; 5.12; 10.24), 
(•) synthetic velocities, (solid line) our theoretical prediction 
(Eq. Q). 



at the end of this letter. Hence, a remains unchanged and 
can be considered as universal. Let us mention that the 
second moment of \ui(k,x)\ is obviously proportional to 
fc" 5 / 3 and (\u £ (k,x)\ 2 ) = (t/L)(\u L (k, x)\ 2 ) for k > 2n/£. 
In the dissipative range (i.e. k > k n ), the flatness appears 
to rapidly increase without any saturation |lCj . This is 
one of the first exp erimental verification of a Kraichnan's 
conjecture [llj,ll2| which is linked to the log- infinitely dis- 
tribution breaking of velocity |l3j and is currently under 
investigations. In the following, we will theoretically es- 
tablish a link between the inertial exponent a and the 
structure function exponent Q : a ~ Ci ~ ^Cz- 

Let us now begin with defining an intermittent (zero 
mean) velocity profile. This is usually done with the help 
of wavelet seri es Il4j . early introduced in the context of 
turbulence [f5lll6l|. 



u(x) 



+oo 2 3 -l 
3 =0 k=0 



^j,k( x ) 



(1) 



where the set {ipj,k(x) = W^^Vx - k)} is an orthonor- 
mal basis of the space of finite energy functions L 2 ([0, L]) 
(see 01) an d ip an admissible "mother" wavelet. The 
wavelet coefficients Cj govern the statistics across scales. 
Generally, the coefficients Cj.fc = 2^^ 2 €j,kPj,k are cho- 
sen as a product of a sign (ej t k — ±1 with equal proba- 
bility) and positive random variables f3j t k that are cho- 
sen so as to be compatible with turbulence longitudi- 
nal velocity statistics, i.e. E(/3? fe ) = 2~ K " (E(.) mean- 
ing here mathematical expectation). Moreover, as al- 
ready predicted by the unifying point of view of Cates & 
Deutsch 0] , statistics of velocity fluctuations are corre- 
lated in space and scale, that can be formalized through 
space-scale correlations of dyadic wavelet coefficients as 

HPtkPVm 



Cfkim ren der additional correlations 



)^JfeTm where the functions 



j,k,l,m 



[|fc2~ J - m2~ l \ + max(2^,2^)] 



(2) 

stating that wavelet coefficients are typically correlated, 
in amplitude, up to the correlation length L. The gener- 
ated velocity profile u cws (x) using Eq. QJ, where (3j t k 
are correlated according to Eq. will be called a 

Cascade Wavelet Series (CWS) 0, HE E3| Using the 
simplest admissible Haar wavelet (i.e. i/j(x) = 1 for 
x G [0;L/2[, ip(x) = -1 for x G [L/2;L[, ^{x) = in- 
stead) and the box for the Short Time Fourier transform 
(i.e. hi(x) = 1 for x G [0;£], and hi(x) = instead), 
it can be shown analytically that E(|2^ ws (fc„, 0)| 2 ) oc 
(£ P /L)k~ 1 ~ < ' 2 , for k„ = Tt2 n /L > ir/£ p and 
L2~ p ((n,p) G N 2 ). Moreover the flatness Fn v (k n ) 
E(|wf ws (fc„,0)| 4 )/(E(|^ ws (fc„,0)| 2 )) 2 behaves as 



/ — 



(1 + C4-2C 2 )(2 + C4-2C 2 ) 



C4-2C2 



(3) 



3 



when l v — > independently of k n . We thus have demon- 
strated that a = (4 — 2(2- Note that the flatness is not 
exactly equal to 2 at the correlation length (£ p = L). 

In order to check our analytical predictions, in par- 
ticular to verify whether our computations depend on 
the box he and the synthesis wavelet ip, we perform a 
statistical study of the process u cws (x), using a more 
regular Daubechies-6 wavelet for the synthesis wavelet 
tp and a Hanning window for he- The method used to 
build the positive random variables 0j h is the classi- 
cal multiplicative cascade model 0, El Ell EU E3 
: recursively, /3 ,o = 1, /3?,2fc = W-^i.fcft-i.fc and 

Pj,2k+i = k Pj-i,k, where the W^\ k (k = I for 

left or r for right) are independent identical ly d istributed 
(i.i.d) positive random variables (see [lil fl7j'). As an 
example, we will study the log-normal case where each 
is a Gaussian random variable of mean //In 2 
and variance a 2 In 2 (leading to the quadratic spectrum 
Cq = i-iq — <J 2 q 2 /2). We have used the set of parameters 
o~ 2 = 0.025 and fi = 1/3 + 3er 2 /2, consistent with ex- 
perimental findings [2]]], so that (2 ~ 2/3, £3 = 1 and 
2C2 ^ C4 = 0.1. Numerically, the infinite sum in Eq. 
is truncated at j = N = 2 18 and is generated over 2 5 in- 
tegral scales. It can be shown that such a stochastic pro- 
cess is not stationary but at first order, its correlation 
function is consistent with Eq. © . We show on Fig. ^b) 
(solid line) the estimation of the flatness of u™ s (k, x) as 
a function of ln(fc//c ma x)- After a k-dependent crossover 
(data not shown) linked, among other reasons, to the ef- 
fect of discretness in the cascade [2(|, the flatness does 
not depend on the wave vector k but significantly de- 
pends on £ p . In Fig. ^c), we have gathered all the 
values of the inertial plateaus using a least-square fit pro- 
viding an error bar estimation. The plateau behaves as 
a power law of the scale l v in accordance with Eq. © 
: C4~2(2 = — 0.1 ± 0.001. The discrepancies between 
the prefactors are mainly linked to the non-stationary 
character of this generated synthetic velocity profile. 

We would like to mention that if wavelet coefficients 
are no longer long-range correlated (take Cfkim = -0 
and if In f3j^ are chosen to be independent Gaussian ran- 
dom variables with mean // In 2-? and variance a 2 In 2-? , 
the corresponding synthetic velocity generated will be 
called a Random Wavelet Series u rws (x) which is inter- 
mittent in a mathematical sense \'2'i . By construction, 
E(|w™ s (fc„,0)| 2 ) = E(\u^ s (k n ,0)\ 2 ), and analytical cal- 
culations performed in the same framework defined in 
the context of CWS (i.e. using a Haar wavelet for tp and 
a box for he) shows that the flatness J- p (k n ) is equal to 
2 for both every wave vectors k n and window lengths 
l p . This property has been checked numerically (with 
the same parameters /1 and <r, and the same synthesis 
wavelet tp and analysis window he as in the CWS case). 
The results are presented in Figs, ^b) and (c). This 
heuristic synthetic process shows that experimental lon- 




FIG. 2: Acoustic scattering experiment in a turbulent jet flow. 

gitudinal velocity data are not only intermittent, but also 
long-range correlated. 

We will turn now to acoustic scattering measurements 
allowing the direct access to a spectral characterization 
of the vorticity distribution. The non-linear interaction 
of an acoustic wave with a turbulent flow gives rise to 
a scattering process of the incident sound wave by the 
turbulent vorticity distribution. As in any scattering ex- 
periment (like e.g. light or neutron scattering), the com- 
plex amplitude of the scattered acoustic pressure field is 
directly related to the incident acoustic amplitude and 
to the spatial Fourier transform of the vorticity distri- 
bution |24j . We have performed such a spectral investi- 
gation on a turbulent axisymetric jet in air, at a Taylor 
based Reynolds number R\ « 300. We use a bistatic 
configuration (Fig. |2J) wherein a plane monochromatic 
acoustic wave with frequency and complex amplitude 
Pinc(t) is directed on the turbulent flow. The complex 
amplitude p SC att{t) of the sound wave scattered at an- 
gle 9 is then recorded along time by a separate acoustic 
receiver. Further details of the experimental apparatus 
can be found elsewhere j^. The scattering process re- 
sults in phase and amplitude modulations of the scat- 
tered acoustic pressure with respect to the incident one: 

Pscatt{t) = (k,tj .pinc(t) where, 

± (k, tj = (ft n ± (x, t) e~ i% - s <fx (4) 
"* ■* Vscatt 

is the spatial Fourier transform of the vorticity compo- 
nent normal to the scattering plane at wave- vector k such 
that \k\ = 47Wo/c x sin {0/2), with c the sound speed. By 
fixing both vq and 9, the scattering experiment allows the 
direct probing, in time, of a well defined spatial Fourier 
mode of the turbulent vorticity distribution character- 
ized by a unique spatial wave-vector k (spectral resolu- 
tion). The price to pay for such a spectral resolution 
lies in some spatial derealization in the physical space 
manifesting itself as a windowed spatial Fourier trans- 
form over a finite volume V sca tt according to equation 
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FIG. 3: Flatness of experimental Fourier modes of vorticity 
vs. ln(£/L), for four inertial wave vectors ki. Dotted line 
slope equal to -0.24. 



The measurement volume is defined by the inter- 
section of the incident and detected acoustic beams and 
mainly depends on 9 and on the size of both acoustic 
transducers. In the present experiment, 9 = 60° and 
the diameter of the circular transducers is 14cm, leading 
to spatial resolutions of order the integral length scale 
of the jet flow L. By varying i/q, at a fixed 9, four 
different wave-vectors ki, i = 1,2,3,4 (in growing or- 
der) have been analysed, corresponding to various length 
scales spanning the whole inertial range of the turbulent 
flow. In the spirit of the first part of this letter, we are 
interested in the influence of the parameter l/L where 
^ ~ (Vscat) 1 ^ 3 is a typical size of the measurement vol- 
ume and L is the integral scale of the flow over the flatness 
^(fci) = (|0^(fc,i)| 4 ) t /(|0^(fc,i)| 2 )|. To this end, we rely 
on the classical selfsimilarity property of the axisymctric 
turbulent jet according to which the integral scale 
L (and all other pertinent scales) increases linearly with 
respect to the distance downstream from the jet nozzle. 
A well known consequence of this statistical selfsimilar- 
ity is the invariability of the Reynolds number for large 
enough distances from the jet nozzle. Several scattering 
experiments have been performed at different distances 
from the jet nozzle, corresponding to different integral 
length scales L. On Fig. [3] is plotted the flatness of the 
modal amplitude Ti(ki) as a function of ln(£/L). Firstly, 
Ti(ki) does not depend on the wave vector ki at a first 
order for I < L, in accordance with Fig. [He). Secondly, 
we see that the flatness behaves as a power law with scale 
I, i.e. F t {ki) ~ P with 7 = -0.24 ± 0.02 when l/L -> 0. 
Analogously with longitudinal velocity, the exponent 7 is 
directly related to the classical exponent of structure 



analysis of Kestener & Arneodo |27| that have been ap- 
plied to a 3D-vorticity field obtained from a Direct Nu- 
merical Simulation (DNS) of Navier-Stokes equations at 
a smaller Reynolds number TZ\ = 140 for which they ob- 
tained C2 - K2 = -°- 22 ± °- 016 - In the same spirit; 
it has been measured that tranverse velocity profiles are 
more intermittent that longitudinal ones [28j. 

To sum up, we have shown that longitudinal velocity 
and vorticity arc intermittent and long-range correlated 
in the physical space thanks to the study of the flatness of 
experimental velocity data and acoustical measurements 
of vorticity Fourier modes. We have seen that vortic- 
ity is much more intermittent than longitudinal velocity. 
Our theoretical study on wavelet series shows that the 
intermittency is intrinsically related to long-range corre- 
lations. We mention the necessity to generalize this ap- 
proach to dissipativc length-scales and wave vectors for 
which log-infinitely divisible principles are violated. 
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experimentale de spectroscopie acoustique multi-echelles) 
and the CNRS. We wish to acknowledge B. Lashermes 
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functions of vorticity considered as a vector field, namely 
7 = — 2£^- It is noticeable that our experimental 
finding is in excellent agreement with a tensorial wavelet 
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